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Abstract. We study the shape parameters of the D-r scalar and vector form factors using 
as input dispersion relations and unitarity for the moments of suitable heavy-light correlators 
i-^ ' evaluated with Operator Product Expansions, including ©(a^) terms in perturbative QCD. For 

the scalar form factor, a low energy theorem and phase information on the unitarity cut are 
^ I implemented to further constrain the shape parameters. We finally determine points on the 

_ real axis and isolate regions in the complex energy plane where zeros of the form factors are 

• excluded. 



1. Introduction 

The knowledge of shape parameters of the Dtt form factors is of interest for the determination of 
the element \Vcd\ of the Cabibbo-Kobayashi-Maskawa (CKM) matrix entering precision tests of 
the Standard Model. The branching fractions of the semileptonic decays D — )• ttIi/ and D — ?• Klu 
have recently been analysed by the CLEO collaboration [H [2] , which renewed the interest in the 
theoretical study of these processes. 

The knowledge on the form factors can be further improved using analyticity and unitarity 
techniques, in particular, the method of unitarity bounds OH]. Using this framework, one can 
obtain bounds on the form factors using as input an integral of the modulus squared of the form 
factor along the unitarity cut, known from unitarity and a dispersion relation for a suitable 
correlator of the same current, which can be evaluated by Operator Product Expansion (OPE) 
in the spacelike region. Using standard mathematical techniques, one can then correlate the 
values of the form factor and its derivatives at different points inside the analyticity domain. A 
review of the method was presented recently in [5]. In this work, we apply the above technique 
to the Dtt form factors. The work is motivated in part by the progress in perturbative QCD 
calculations, which now yield the heavy-light correlators of interest to order [6] for a massless 
light quark and partly form the new experimental analyses of the semileptonic D decays. 

We investigate the shape parameters entering the Taylor expansion around t = 0, 

fkit) = h{0) (l + X',^ + ^A^'^ + . . . ^ , A; = +, 0, (1) 



and derive allowed ranges for the slopes and the curvatures A^. Here, and fo{t) are 

respectively the vector and scalar form factors and are related to each other by the relation, 

fo{t)=U{t)+ ^,J_^^J .{t), t = q^ = {p-p'f. (2) 

We further find points on the real t-axis and isolate regions in the complex t-plane where 
zeros of the form factors are excluded. The knowledge of the zeros are of interest, for instance, 
for the dispersive methods based on phase (Omnes-type representations) and for testing specific 
models of the form factors. 



2. Reviev^ of the method 



To start with, we have the heavy-light invariant amplitudes n_|_(g^) and IlQ{q'^) defined by the 
vector-vector correlation function, 

- {q^g^" - q^q^Ii+iq^) + gVno(g^) = i j <fxe'^''{Q\TV^'{x)V''H<d)\Q), (3) 

where = d'y^c. 

We then consider the moments of the invariant amplitudes at tf' = 0, 

^^"'";!!(^[^^(^'^].-o' ^ = +'0' (4) 

which satisfy dispersion relations of the form, 

(„) 1 p Imnfc(t+^ 
Xfc =-/ dt , A; = +,0, (5) 



vr 
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where t± = [Md ± M^f . 

It is well-known from QCD that the amplitude n_|_(g^) satisfies a once subtracted dispersion 
relation, while for no((7^) an unsubtracted relation converges. The quantities and Xo^ are 
hence defined for n > 1 and n > 0, respectively. These relations are connected to the form 
factors and /o(t) by means of unitarity, 

lmn„(* + .)>^^l(i^^±>^^^|/„(*)|^ (7) 

respectively for the vector and scalar form factor, which hold for t > t+. The above expressions 
are the unitarity sums for the spectral functions, restricted to the contribution of the Dn states 
in the isospin limit. 

(n) 

OPE on the other hand allows for the calculation of Xk ^.s the sum of perturbative and 
nonperturbative contributions. While the perturbative parts of the moments of heavy-light 
correlators for n < 7 were calculated up to order in [B], the leading non-perturbative 
contribution of the quark and gluon condensates were obtained from [7], the corresponding 
expressions of which are presented in [8]. For the vector and scalar form factors, we specifically 

use the dispersion relations for the moments x+\x+\x+^ and Xo'\ x'o\ X(j'^ respectively, as a 
result of which we obtain, for each form factor, a family of three different constraints. The final 



allowed domain for the parameters of interest will be the intersection of the three individual 
domains. 

Prom the dispersion relations ([5]) and the unitarity conditions given in eqns.(l6]) and ([7]), it 
follows that each form factor fk{t) satisfies a set of integral inequalities written as, 

-/ dtp^^\t)\Mt)f<xP, k = +,0, (8) 

where the weights p^^\t) are the product of with the phase space factors entering the 

unitarity relations. Using the fact that the form factors and /o(0 are analytic functions 

in the complex t-plane cut along the real axis from 4+ to oo, we apply standard techniques to 
derive from eqn.([8]) constraints on their values, in particular on the shape parameters and on the 
regions in the real and complex energy plane where zeros are excluded. Specifically, we cast the 
problem into a canonical form via a conformal map, and derive a determinant which is central 
to our investigations for obtaining bounds on the shape-parameters O [HI [9] and for finding the 
exclusion regions of the zeros. The quantity /fc(0) defined in eqn.([T]) goes as an important input 
to our work. Besides, additional information for the scalar form factor provided by a low-energy 
soft-pion theorem namely the Callan-Treiman relation [TU] , is exploited to obtain more stringent 
bounds on the shape parameters of the scalar form factor For the Dtt case, the corresponding 
expression is given below, 

fo{Ann) = fo/fn, (9) 

where Ad^^ = M|, — is the relevant Callan-Treiman point and and are the meson 
decay constants. 

Even more stringent constraints on the shape parameters can be obtained if we also have 
some information on fk{t) on the unitarity cut, in particular, if the phase 6k{t) defined as, 

fkit + ie) = \fkit)\e''''^'\ k = +,0, (10) 

is known along a low-energy interval, t+ < i < tin- This information on the phase can be 
implemented using the technique of generalized Lagrange multipliers and involves the solution 
of an integral equation. A review of the method and more references are given in [5]. 



3. Inputs and Results 



The perturbative and nonperturbative contributions to the moments Xk have been tabulated in 
[S]. The value of the inputs used for /^(O) and fo{Ann) which reads 0.67+[];J? and 1.58 ± 0.07 
respectively are obtained from a recent work based on analyticity [TT]. We roughly estimate the 
phase of the Dtt form factors from the masses and widths of the resonances dominant at low 
energies, namely, from the relativistic Breit-Wigner parametrization given by. 



6(t) = arctan ( — ^ 

\Ml-t) 

where Mn is the mass and T[t) the energy-dependent width 



i'4i^). (") 



written in terms of the angular momentum J, the width and the centre of mass momentum 
q{t) = y^{t-t_){t-t+)/4t. 



The central values of the masses and widths of the lowest charged Dn vector and scalar 
resonances, namely the D* and Dq respectively are listed in [12]. We note that the vector 
resonance D* is very close to the threshold so that a reasonable expression of the phase cannot 
be obtained from a Breit-Wigner parametrization, but in the scalar case we can assume that the 
phase (^o(^) is reliably described by the expression ()lip with J = 0. For our analysis, we consider 
the phase up to the point tin = (2.6 GeV)^, close to the first inelastic Dr] channel opening at 
2.42 GeV. 
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Figure 1. Constraints on the slope and curvature of the vector form factor obtained using as 
input different moments of the correlator. The point indicates the slope and curvature of the 
pole ansatz given in [13j . 





Figure 2. Constraints on the slope and curvature of the scalar form factor obtained with the 
moment Xo \ from the standard bound and by including the phase and the low-energy theorem. 
The point indicates the slope and curvature of the pole ansatz given in [13j . 

The results of our analysis are presented in FiglH where the the interior of the ellipses 
represent the allowed domains in the slope-curvature plane for the vector form factor, obtained 
with three moments of the vector correlator with just /4-(0) and no other input. While the 
smallest or in other words the best domain is obtained from the lowest moment, we find that 
the higher moments tend to slightly reduce the domain, since one must take the intersection of 
all the domains in order to fulfill simultaneously the constraints. We also indicate the slope and 
curvature of the simple pole ansatz [T3], which is seen to satisfy the unitarity constraints coming 

from our analysis. In FigjJl we show the bounds derived from the lowest moment Xo''^ of the 
scalar form factor where we also use information on the phase and the soft pion theorem. The 
effect of these additional constraints can clearly be seen in the figure. We also show the slope 
and curvature from the pole ansatz [13] which satisfies our constraints. It may be noted that 



the above bounds were obtained with the central values of the input parameters. By varying 
simultaneouly all the input values we can obtain more conservative regions, which are slightly 
larger than the domains shown in the figures. 
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Figure 3. Domain without zeros for the vector form factor, obtained from the lowest moment 
and the input /+(0). 

We now present results on the zeros of the scalar and vector Dn form factors. We have 
considered in the following only the constraints imposed by the lowest moments, corresponding 
to n = 1 in the vector case and n = in the scalar one, using /+(0) as input quantity. For the 
vector form factor, real zeros are excluded in the range (-1.0, 0.80) GeV^ while for the scalar 
case it is (—2.51, 1.55) GeV^ for the standard bounds. With the inclusion of the low-energy 
constraint in the case of the scalar form factor, the range is increased to (-3.55, 3.89) GeV^. We 
present in Figs. [3] and S] the regions of excluded zeros in the complex t-plane for the vector and 
scalar form factors respectively. In the scalar case the zeros are excluded in a larger domain if 
the low-energy constraint is also imposed. 



E 



■10 



Figure 4. Domain without zeros for the scalar form factor, obtained from the lowest moment 
Xq^'' and the input /o(0) (smaller region) and using in addition the CT constraint (larger region). 



4. Conclusions 

In this work, we have used analyticity and unitarity techniques to derive information on the 
Dtt form factors, which might be of interest for the determination of the element \Vcd\ of the 
CKM matrix. We have applied the formalism of unitarity bounds and have derived a family of 
constraints to be satisfied by the shape parameters of both the scalar and vector form factors at 



t = and explored regions on the real axis as well as in the complex plane where the form factors 
cannot have zeros. The exclusion regions for the zeros that we have isolated basically cover a 
significant part of the entire low energy region. They are of practical interest, for instance, for the 
dispersive methods based on phase (Omnes-type representations) and for testing specific models 
of the form factors. The technique can be extended also for deriving new parameterizations 
for the Dtt form factors in the semileptonic region, similar to those proposed in |14j for the 
semileptonic B ^ ttIv decay. The new parametrizations properly implement the singularities 
related to the lowest charmed resonances and are useful for the description of the semileptonic 
data (see Ref.[5] for more details). 
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